Self diffusion of particles in complex fluids: temporary cages and permanent barriers 
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We study the self diffusion of individual particles in dense (non-) uniform complex fluids within 
dynamic density functional theory and explicitly account for their coupling to the temporally fluc- 
tuating background particles. Applying the formalism to rod-like particles in uniaxial nematic and 
smectic liquid crystals, we find correlated diffusion in different directions: The temporary cage 
formed by the neighboring particles competes with permanent barriers in periodic inhomogeneous 
systems such as the lamellar smectic state and delays self diffusion of particles even in uniform 
systems. We compare our theory with recent experimental data on the self diffusion of fluorescently 
labelled filamentous virus particles in aqueous dispersions in the smectic phase and find qualita- 
tive agreement. This demonstrates the importance of explicitly dealing with the time-dependent 
self-consistent molecular field that every particle experiences. 

PACS numbers: 66.10.cg, 47.57.J-, 61.20.Gy, 61.20.Lc 



Phenomena such as multi-site hopping in microstruc- 
tures, void diffusion in colloidal crystals, dynamical het- 
erogeneities of colloidal glasses and the self-assembly of 
micelles, supramolecular polymers and viruses are strik- 
ing examples of the intriguing dynamics of complex flu- 
ids, whose understanding remains relatively rudimentary 
despite intense research spanning many decades. The 
main reason for this state of affairs presumably is an in- 
complete understanding of temporal fluctuations of the 
fluid structure and how these couple to single-particle dif- 
fusive motion. To circumvent this problem, self diffusion 
of particles in locally or globally inhomogeneous fluids 
is therefore often studied in a fixed background poten- 
tial mimicking the actual fluid structure but completely 
ignoring the inherently fluctuating nature of it [l| . 

It is the aim of this Letter to theoretically investigate 
the influence of the local fluid structure and the tem- 
porally fluctuating background on the self diffusion in 
uniform and non-uniform complex fluids. As a simple 
yet quite interesting example we apply our treatment, 
based on dynamical density functional theory, to disper- 
sions of elongated colloidal particles in uniaxial nematic 
and smectic liquid-crystalline phases. In the former the 
orientational degrees of freedom are frozen out and in the 
latter also one positional degree of freedom. This allows 
us to compare our theory with results of a very recent ex- 
perimental study of the unusual self diffusion in aqueous 
dispersions of the filamentous bacteriophage fd Q . 

As we shall see, the local fluid structure forms a tempo- 
rary cage around every test particle that initially hinders 
its free self diffusion but that decays at later times. Re- 
markably, this temporal caging effect of the background 
particles can produce a coupling between motion in dif- 
ferent directions in particular if the fluid is symmetry bro- 



ken. Both these phenomena influence the self diffusion 
in structured fluids and cannot be accounted for by pre- 
suming a fixed molecular background field. In fact, even 
if a fluctuating molecular fleld was presumed, its effect 
can only be predicted beyond the usual linear analysis of 
density fluctuations, i.e., they are inherently non-linear. 

Focal point of our discussion are the Van Hove corre- 
lation functions that probe diffusive processes [3|. We 
generalize the formalism introduced for simple fluids in 
Ref. [j|. The key idea is to deflne conditional densities 
for which equations of motion can be prescribed [J] and 
which for non-uniform or complex fluids are non-trivially 
related to the well-known Van Hove correlation functions. 
The formalism, although applied here to athermal sys- 
tems can be generalized straightforwardly to any sys- 
tem for which a free energy functional can be written 
down, including, say, thermotropic liquid crystals within 
Landau-de Gennes theory 5j. 

Here, we only outline the main ingredients of the the- 
ory. Consider an equilibrium fluid of N particles with 
arbitrary degrees of freedom x = (r,uj, . . .) and define 
self (s) and distinct (d) conditional densities as 

a(x,i|x',0):=-^^^<5(x-X„(i))<5(x'-X„(0))^,(l) 

1 / ^ 
Cd(x,t|x',0):=— -( J2 '5(x-X„(t))<5(x'-X„,(0)) 

n,n —1 

respectively, where (• • • ) is the average over all equilib- 
rium trajectories, 6 the Dirac-i5 in configuration space, 
X„(t) the configuration of particle n at time t, and £i(x') 
the one-particle equilibrium density at configuration x'. 



Note that the distribution of equiUbrium trajectories is 
time-translationally invariant. The conditional densities 
Cg.d and the Van Hove self and distinct correlation func- 
tions Gs,d \m ^^'^ related via 

G„d(Ax, f) = ^ Jdx' Qd(x' + Ax, t|x', O)g(x'), (2) 

where x' + Ax denotes the displacement of configuration 
x' by a suitably defined offset Ax in configuration space. 
According to Eq. ^, the conditional densities at zero 
time t = read Cs(x,0|x',0) = (5(x - x') =: ni(x,0) 
and Cd(x, 0|x',0) = p(x)f/(x^x') ~: nu(x, 0) with g the 
pair distribution function [3| and ni.u the one-particle 
equilibrium densities of a fluid of one labeled (1) and 
N — 1 unlabeled (u) particles in which the labeled par- 
ticle is fixed in configuration x'. Upon releasing the 
fixed labeled test particle at time t = 0, the neighbor- 
ing host fluid relaxes towards a new equilibrium state 
and it is assumed that the time-dependent two-particle 
correlators Cs,d(x, t|x', 0) equal the one-particle densities 
ni^u(x, i) for all times t > 0. This identification is remi- 
niscent of Onsager's regression hypothesis and implies the 
fluctuation-dissipation theorem to hold, which is strictly 
proven only within linear response theory ^] . A free en- 
ergy density functional F[ni,nu] [7| describing the fluid 
automatically produces the one-particle equilibrium den- 
sity distribution g and the pair distribution function g as 
well as expressions for the local chemical potentials 



//i,u(x,t, [ni,nu]) := 



5F 



Sni.uix) 



(3) 



ni(- ,t) ,n^^{• ,t) 



that within a generalized Fickian approximation give rise 
to the equations of motion 8] 



9 , . 

^?lKu(x,i) 



(4) 



Vx • (n-Lu(x, t)r(x, i, [ni, nj) • VxMKu(x, i, [ni, nj)) 

with r the mobility matrix. Note that the Fickian ap- 
proximation leading to Eq. (|4]) implies the hypothesis of 
local equilibrium and the neglect of hydrodynamic inter- 
actions, i.e., the free draining limit to effectively hold. 
In dense systems the latter presumption seems reason- 
able due to screening of hydrodynamics. The mobil- 
ities are then to be interpreted as renormalized ones. 
Integrating this closed set of equations for m^u sub- 
ject to the specified initial conditions and interpreting 
ni^u as conditional densities Cg.d, one obtains the Van 
Hove correlation functions Gs,d from Eq. Q. For a 
uniform, simple fluid, where the particles possess only 
translational degrees of freedom, x = r, and the one- 
particle density g is spatially constant, the conditional 
densities C's,d(r, t|r', 0) depend only on r — r', so that 
Gs,d(Ar,i) = Gs,d(Ar,i|O,0) = ni,u(Ar,t) and Eqs. © 
and (HI can be interpreted as equations of motion for the 



Van Hove correlation functions Gs,d and avoids having to 
explicitly deal with kinetic equations for conditional den- 
sities that in actual fact are two-point correlators. Note 
that the quantities called "van Hove functions" in Ref. 4| 
are actually conditional densities, which, in general, dif- 
fer from the well-known Van Hove correlation functions. 
The general conditional density formalism described here 
is applicable to any non-uniform or complex fluid and 
offers a route to analyze the relaxational dynamics of a 
wide range of interesting complex fluids. 

Motivated by the very recent measurements of the Van 
Hove self correlation function of aqueous solutions of the 
bacteriophage /d [2] , a filamentous virus particle of about 
900 nm length and 7 nm width, we consider a free energy 
functional F[ni,nu] describing a lyotropic liquid crystal 
of (stiff) hard rods of length L and diameter D with 
L ^ D. As the root-mean-squared angle between the 
axis of a rod and the director scales as 0{D/L) in this 
limit [9| , the orientational degrees of freedom can be ig- 
nored, so the rods are assumed to be oriented parallel 
to the director, which itself is assumed to be fixed. This 
means that the model particles possess only translational 
and no orientational or internal conformational degrees of 
freedom such as arising from a bending flexibility. For the 
free energy density F[ni, riu] we for reasons of simplicity 
invoke the second virial approximation 'lO]. Although 
not accurate it is known to capture the main features 
of the structure of the smectic phase near the nematic 
transition point, which suffices for our purposes. Our 
model of perfectly parallel hard rods is simple but not 
overly simple, as quantitative precision is not required 
and the neglect of features such as particle flexibility, 
higher virial terms, or (screened) hydrodynamic interac- 
tion can in principle be accounted for by renormalization 
of the model parameters. The bare translational mobil- 
ities Fii and r_L parallel and perpendicular to the direc- 
tor give rise to the parallel translational diffusion time 
T :— /JL^/Fii and the diffusion rate ratio 7 := 't/lT- 

We calculated the Van Hove correlation functions 
Gs,d{z,r,t) as a function of the parallel displacement z, 
the perpendicular displacement r, and the time t. Re- 
sults for the largely arbitrary but definitely representa- 
tive choice of parameters 7=1 (representing fast radial 
diffusion) and 7 = lO^"' (slow radial diffusion) for the 
three cases of a nematic state N at a chemical potential 
/3A/i = —0.388 relative to that at the nematic-smectic 
transition, a weakly smectic state Si with /3A/i — 0.612, 
i.e., at a density just above the smectic transition, and 
a strongly smectic state S2 with /3A/i = 2.612 are shown 
in Figs. [11 [2l and [3l In passing we note that paral- 
lel diffusion becomes independent of the parameter 7 if 
7 > 1, because perpendicular relaxation is then much 
faster than parallel relaxation according to our numeri- 
cal results. The interlayer distances of the smectic phases 
within the present model are 1.4L for Si and 1.3L for 




FIG. 1: Radially integrated Van Hove self correlation function 
Ga{z,t) as a function of the axial displacement z in units of 
the rod length L and time t in units of the diffusion time r for 
diffusion rate ratio j — 1 for (a) the nematic state N, (b) the 
weakly smectic state Si, and (c) the strongly smectic state 
S2. The thin vertical lines in (b) and (c) indicate the centers 
of the smectic layers. 



^2; these values would be somevsrhat smaller had we in- 
cluded particle flexibility 11 1 and higher virial terms [ld| ■ 
For these two smectic states, the calculated lamellar den- 
sity undulations correspond to a self-consistent molecular 
field barrier height of 2.1kBT and A-lksT, respectively. 

The radially integrated Van Hove self correlation func- 
tions Gs{z,t) := 2Tr J^drrGs{z,r,t) presented in Fig. [1] 
exhibit a spatial broadening due to self diffusion, which is 
slower the larger the background potential barrier height 
is. Whereas Gs{z,t) is a concave function of z for the 
nematic state N for all times t, shoulder peaks develop 
in the smectic states Si and S2 located at the centers of 
the smectic layers in accord with the experimental find- 
ings of Ref. [4] . These peaks are the manifestation of the 
existence of the average self-consistent field due to the 
equilibrium one-particle density. The curves of the ne- 
matic state and the envelopes to the curves of the smectic 
states are not Gaussians (see Fig.[3}3) because of the influ- 
ence of the background fluid that cages the test particle. 
Upon approaching the nematic-smectic phase transition 
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FIG. 2: Cuts of the Van Hove distinct correlation function 
Gdiz,r,t) as a function of the axial displacement z in units 
of the rod length L, the radial displacement r in units of the 
rod diameter D, and time t in units of the diffusion time t 
for diffusion rate ratio 7 = 1 for (a,c) the nematic state N 
and (b,d) the weakly smectic state Si. The corresponding 
curves for the strongly smectic state S2 are similar to those 
in (b,d). The thin vertical lines in (b) indicate the centers of 
the smectic layers. 



from the high-density, smectic side we found an increase 
of the time that it takes for the first shoulder peak to 
appear (see also the curves t — 0.3r in Figs. [T]d and c). 
This we attribute to influence of the critical slowing down 
of the collective dynamics of the host fluid when nearing 
the smectic spinodal. This may well explain the absence 
of shoulder peaks in Fig. 2b of Ref. y presuming that 
this time exceeds the time of measurement. 

Figure [2] shows the time-evolution of the structure of 
the fluid background for the nematic and the smectic 
phase in terms of the Van Hove distinct correlation func- 
tion Gd{z, r, t). The initial Van Hove distinct correlation 
function Gd{z,r,t — 0) coincides with the pair distribu- 
tion function and the nearest-neighbor (solvation) shell 
in the vicinity of the origin corresponds to a cage around 
the labeled test particle. This cage initially slows down 
the self diffusion of the labeled test particle, but it dis- 
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FIG. 3: (a) Parallel mean-squared displacement (z^) in units 
of the squared rod length I? and (b) Binder cumulant U 
(see main text) of the labeled particle as a function of time 
t in units of the diffusion time r for the diffusion rate ratios 
7 = 1 and 7 = 10^* for the nematic state N, the weakly 
smectic state Si and the strongly smectic state S2. The points 
represent the results obtained within the present fluid model, 
whereas the lines describe the one-dimensional self diffusion 
in a temporally fixed background potential corresponding to 
the equilibrium one-particle density. 



solves with time. Eventually, the labeled particle diffuses 
effectively only in the permanent background of the unla- 
beled particles. The presence of the temporary cage that 
adds to the permanent equilibrium background barriers 
between the particle layers is a direct consequence of the 
local fluid structure, which cannot be accounted for by 
fixed background models as in Ref. [l|]. For 7=1, the 
relaxation of G^iz^ r, i), i.e., the decay of the temporary 
cage, takes place on a time scale of the diffusion time 
T. If the perpendicular diffusion is strongly suppressed, 
e.g., for 7 — 10^^, this relaxation takes place on a larger 
time scale, indicating that the dissolution of the tempo- 
rary cage is linked to the perpendicular diffusion. Hence, 
the time-dependent fluid background mediates a coupling 
between parallel and perpendicular self diffusion. 

The influence of the fluid structure on the self dif- 
fusion of individual test particles can be inferred from 
Fig. [3] that shows our results for their mean-squared 
displacement (z^) along the director of the fluids and 
for the Binder cumulant U := (z )/(3(2;^)^) — 1, which 
quantifles deviations from Gaussian behavior []|, where 
(z*^) :— Jdz z'^Gsizjt). The significant difference be- 
tween results obtained within fixed background models 
[l| and our fluid background model highlight the impor- 
tance of the temporal nature of fluid structure in gen- 
eral and of the sideways diffusion in the intermediate- 
time regime of uniaxial fluids in particular. For each 



of the three states investigated, one spatially homoge- 
neous and two spatially inhomogeneous, caging causes 
diffusion initially to be slowed down giving rise to a non- 
vanishing Binder cumulant. The crossover from early to 
late-stage diffusive motion of test particles in a congested 
fluid causes this deviation from Gaussian behavior and 
can be linked to potential barriers that can be tempo- 
rary or flxed. In either case, they are due to correlations 
between particles in the fluid that cannot reasonably be 
described by a flxed potential. Our calculations show 
that in the crossover region sub-diffusive behavior may 
be inferred over a limited time range, as done in Ref. [2]- 
The slowing down of fd virus particle diffusion in the 
smectic phase was rationalized in Ref. [2| by considera- 
tions involving a fixed periodic molecular field only. In 
our view this does not do justice to the complexity of the 
problem that involves a coupling of between-layer and 
in-layer diffusive processes. This can be inferred from a 
comparison of the cases 7 = 1 and 7 — 10~^ for the states 
N and Si where a reduction of the perpendicular diffu- 
sivity gives rise to a slower parallel self diffusion. On the 
other hand, the mean-squared displacement for strongly 
smectic states for which the permanent potential barriers 
are quite large is virtually independent of 7, as the initial 
delay of the diffusion is due to the high, 7-independent 
permanent barriers and not due to the temporal cage. 

In summary, we found a remarkable influence of the 
fluctuating local structure on the self diffusion of parti- 
cles in spatially (non-)uniforni complex fluids. Diffusion 
delays initially even in uniform fluids due to a tempo- 
rary cage formed by neighboring particles. This coupling 
of the motion of a test particle to its surrounding parti- 
cles causes motion in different directions to become cou- 
pled, in particular in anisotropic fluids. Fixed molecular 
background models cannot describe this and have to be 
replaced by more sophisticated approaches such as the 
dynamic density functional formalism presented here. 

We thank M. P. Lettinga, E. Grelet, and B. Vorse- 
laars for useful comments as well as M. Schmidt, A. J. 
Archer, and P. Hopkins for a clarifying communication. 
This work is part of the research program of the "Sticht- 
ing voor Fundamenteel Onderzoek der Materie (FOM)" , 
which is financially supported by the "Nederlandse Or- 
ganisatic voor Wetenschappelijk Onderzoek (NWO)" . 



* Electronic address: 'm.bier@uu.nl 
[1] H. Risken, The Fokker-Planck Equation (Springer, 

Berlin, 1996); B. Vorselaars, A. V. Lyulin, K. Karatasos, 

and M. A. J. Michels, Phys. Rev. E 75, 011504 (2007). 
[2] M. P. Lettinga and E. Grelet, Phys. Rev. Lett. 99, 197802 

(2007). 
[3] J.-P. Hansen and McDonald, Theory of simple liquids 

(Academic Press, London, 1986). 
[4] A. J. Archer, P. Hopkins, and M. Schmidt, Phys. Rev. E 



75, 040501 (2007). [8] W. Dieterich, H. L. Frisch, and A. Majhofer, Z. Phys. B 
[5] P.-G. de Gennes and J. Prost, The physics of liquid cry- 78, 317 (1990). 

tals (Oxford University Press, Oxford, 1995). [9] T. Odijk, Macromoleculos 19, 2313 (1986). 

[6] U. Marini Bettolo Marconi, A. Puglisi, L. Rondoni, and [10] B. Mulder, Phys. Rev. A 35, 3095 (1987). 

A. Vulpiani, Phys. Rep. 461, 111 (2008). [11] P. van der Schoot, J. Phys. II France 6, 1557 (1996). 
[7] R. Evans, Adv. Phys. 28, 143 (1979). 



